Biochimica et Biophysica Acta (BBA) - Biomembranes

Lipid packing defects expose hydrophobic regions of membranes and influence protein-membrane recognition, membrane remodeling, and lipid sensing. Quantitative analysis of these defects remains challenging for heterogeneous membranes with strong curvature, mixed bilayer and monolayer organization, and substantial neutral-lipid content. Here, we present Lipid Defect Analyzer (LipiDeA), an open-source Python framework for quantitative analysis of lipid packing defects in complex membrane systems. LipiDeA reconstructs membrane surfaces using triangulated meshes that adapt to local membrane geometry and identifies packing defects from lipid coverage of the reconstructed surface. The framework further classifies bilayer and monolayer regions, resolves defect properties according to local membrane curvature, and distinguishes defects associated with phospholipid tails and neutral lipids. Validation using planar and buckled DPPC bilayers reproduced previously reported packing defect constants and showed that defects preferentially occur in positively curved membrane regions. Application to lipid-droplet-like membranes revealed composition- and topology-dependent defect formation, with the largest defects occurring during the transition from bilayer-dominated to monolayer-containing membranes and distinct defect properties for DAG- and TAG-containing systems. These results highlight the importance of local membrane geometry and composition in determining lipid packing defects.

Amyloid {beta}-protein (A{beta}), which is known for its role in Alzheimer's disease (AD), is also involved in several physiological functions. Of the two predominant isoforms, A{beta}40 and A{beta}42, the latter is more aggregation-prone and more active in health and disease. In the brain, A{beta} interacts with free lipids, which may affect its assembly pathways. Here, thioflavin T (ThT) fluorescence spectroscopy is combined with molecular dynamics (MD) to investigate A{beta}40 and A{beta}42 co-assembly with free dimyristoylphosphatidylcholine (DMPC) and DMPC/cholesterol (CHOL) lipids. ThT measurements at several protein-to-lipid (PL) ratios show isoform- and lipid composition-dependent kinetics of A{beta} fibril formation with clear isoform-dependent lipid effects at a PL ratio of 1:5. Early co-assembly stages of six A{beta}40/A{beta}42 with DMPC and DMPC/CHOL lipids at this PL ratio are then examined by MD. MD results indicate that lipids alter the secondary structure more and interact more avidly with A{beta}40 than with A{beta}42. In turn, A{beta}40 increases the planar order of lipid bicelles more than A{beta}42. Unlike A{beta}40, A{beta}42 decreases the nematic order of CHOL. In DMPC/CHOL lipids, A{beta}42 but not A{beta}40 exhibits two A{beta}-lipid interaction modes associated with mutually exclusive A{beta}42-lipid and A{beta}42-A{beta}42 interactions. Our MD findings provide a plausible explanation for an experimentally observed increase in the lag time for A{beta}40 fibril formation in DMPC lipids and even more so in DMPC/CHOL lipids that contrasts with much smaller effects of lipids on A{beta}42 fibril formation and offer additional atomistic-level insights into A{beta} co-assembly with lipids amenable to experimental testing.

Small multifunctional viral membrane proteins called viroporins are attractive drug targets but are difficult to structurally characterize in a lipid-bilayer environment. We have successfully reconstituted the full-length Envelope (E) protein from SARS-CoV-2 in a lipid-bilayer, co-assembled with a native host factor that binds to E, the PDZ domain 2 of ZO1. We have leveraged the enhancement afforded by dynamic nuclear polarization NMR to determine that the ensemble-averaged oligomer state of E is a pentamer. Complexing the protein with a host PDZ domain did not impact the oligomer state, but resulted in structuring of the last four residues of E, and immobilization of the PDZ domain near the lipid surface. Combined with secondary structure information derived from the chemical shifts, these data were used to model an E protein pentamer together with the PDZ protein dimer. Inclusion of a relevant host factor, as well as a complex lipid bilayer that mimics conditions in the ER-Golgi intermediate compartment, approaches native conditions, and addresses the ongoing controversy regarding the functional oligomer state in an infected cell.

Lipid nanoparticles (LNPs) are a leading drug delivery platform for mRNA or other nucleic acid-based therapeutics. LNP platforms are typically composed of ionizable aminolipids, phospholipids, sterols, and lipids tagged with polyethylene glycol. The ionizable aminolipids are typically cationic at low pH, allowing them to encapsulate negatively charged mRNA. The pH at which the ionizable head groups change charge states, the apparent pKa, is crucial to understand the characteristics and formulation of LNP. However, this is difficult to measure as the apparent pKa differs from the pKa of the aminolipids in bulk solution due to the complex chemical and charge environment. In this study, we use molecular dynamics (MD) simulations with alchemical free energy perturbation to calculate the pKa shift as a function of lipid composition, treating LNP-like and vesicle-like environments. Our results show that the pH-dependent structure of the systems affects the trends in pKa shifts in complex ways, driven by factors including surface electrostatics and lipid internalization. Our all-atom simulation protocol allows us to efficiently estimate the apparent pKa of LNP-like models formed with novel ionizable aminolipids and estimate how it will change with composition.

Liquid-liquid phase separation (LLPS) of biomacromolecules has been shown to influence a large variety of biological processes via the formation of membraneless organelles. Aberrant LLPS, however, has been implicated in a number of oncological and neurodegenerative disorders. Proteins related to Alzheimer's, Parkinson's diseases and other major amyloidoses have all been shown to form these condensates. Interestingly, some amyloidogenic proteins were demonstrated to form heterotypic droplets, where both molecules coexisted within the same condensate and affected each other's aggregation. While such cross-interactions have been reported numerous times under non-LLPS conditions, the heterotypic droplet assemblies were only shown for a very limited number of amyloidogenic proteins. In this work, we examine heterotypic LLPS of prion proteins and the pro-inflammatory S100A9, whose cross-interaction under non-LLPS conditions we confirmed in a previous work. The results show that this interaction under condensation-promoting conditions enhances their aggregation, leading to the formation of a higher concentration of amyloid structures. These findings suggest that the pro-inflammatory S100A9 may have a critical role in promoting prion protein aggregation via the assembly of heterotypic condensates under high molecular crowding conditions.

Carbohydrates are essential molecules involved in a wide range of biological processes. While their structural representation in the Protein Data Bank (PDB) has been hindered due to their highly diverse and flexible nature, advancements on X-ray crystallography and CryoEM are accelerating the obtaining of their experimental structures, which in turn supports improved analyses of their structural properties, such as the role of their non-chair conformations on biological processes. In this context, the current work provides a comprehensive and up to date overview of carbohydrate structures within the PDB through classification, quantification, and characterization of monosaccharides and glycan chains, incorporating Cremer-Pople puckering coordinates to evaluate pyranose conformations. Our findings show that, even with a steep increase in CryoEM structures including more complex oligosaccharides, carbohydrate data remains concentrate dinsmallglycanchains, with a about 4% of non-chair conformation sintopresolution structures. Free energy landscapes for the pyranose ring in the presence and absence of nine different binding proteins, including both enzymes and non-enzymes, reinforce the preferred stability of chairs and the rarity of conformational induction upon binding.

Finite-bandgap charge-based FET biosensors are known to be insensitive to electronic tunneling noise in the semiconductor and interference from weakly-charged proteins in physiological samples. However, their signal is still corrupted by pH fluctuations and condensing counterions in the bulk electrolyte. Here, we demonstrate that electrolyte noise is screened by the Stern layer beyond the linear Debye-Huckel limit, when the interfacial potential drop is larger than the thermal voltage. Using a fully nonlinear Gouy-Chapman-Stern (GCS) theory, a universal calibration curve is derived and shown to collapse the low-noise FET signals for miRNA, extracellular vesicles (EV), and protein, captured either directly or through silica and magnetic nanoparticle charge reporters, over eleven decades of analyte concentration. The theory involves the dissociation constant KD and a gain parameter {chi}, which represents the competing effects of the reporter charge and buffer Debye screening that must be large to produce the low-noise GCS condition. The limit of detection (LOD) is shown to be KD/{chi}, which can be optimized by varying the probe density, charge reporter, and buffer solution. An optimized hybridization buffer then yields {chi} ~ 22 (LOD ~1 pM) for miRNA detection, and multivalent capture reduces KD by four orders of magnitude for EV detection (relative to protein detection) with {chi} ~ 34 (LOD ~10 fM).

Dictyostelium discoideum is a genetically tractable amoebozoan model organism widely used to study conserved eukaryotic cellular processes, yet a comprehensive map of its subcellular proteome has been lacking. Here, we combined subcellular fractionation with label-free mass spectrometry to resolve a global protein localization map for D. discoideum. We detected 6,337 proteins (~50% of the predicted proteome) and generated protein abundance profiles across 10 subcellular fractions. Using a curated marker set and a support vector machine classifier with a median cutoff, we assigned high-confidence localizations to 3,169 proteins across 17 subcellular compartments. Independent validation using sequence-based targeting predictions strongly supported the compartment assignments. Together, this study provides the first global subcellular localization map of the Dictyostelium proteome, establishing a foundational resource for functional, cell biological, and evolutionary analyses. We are currently preparing manuscripts that further analyze (i) the membrane trafficking system of Dictyostelium, including detailed investigation of the contractile vacuole, Jotnarlogs, and patchy proteins, and (ii) the composition of the Dictyostelium peroxisome, with a focus on the subcellular localization of sterol biosynthesis enzymes. Investigators interested in using these data are encouraged to contact us prior to publication to avoid overlap and to facilitate coordinated and collaborative use of this resource.

G protein-coupled receptors (GPCRs) represent the largest membrane protein family in the human body and the major targets of FDA-approved medications. Unlike some other physiological proteins, GPCRs are highly plastic and often exist as collections of conformational ensembles in the environment to exert their function. Therefore, understanding how a GPCR adopts different conformational ensembles in response to ligands and signaling partners is essential for drug development. 19F-NMR has emerged as one of the most promising tools in this regard; however, a well-established 19F-labeling system is still lacking, including an effective 19F probe that meets both environmental susceptibility and NMR sensitivity requirements, even though progress has been made dramatically in recent years. Thus, this manuscript reports a novel conformational probe (4-(bromomethyl)-ethoxy-2-(trifluoromethyl)benzene, BTMB) that offers paramount benefits for future GPCR conformational dynamics studies. To rationalize the mechanism behind BTMB's superiority for future probe development, we performed MD simulations, referenced to the current state-of-the-art probe, BTFMA. Our study suggests that the spatial proximity of aromatic residues significantly influences the chemical shift dispersion and linewidths of related 19F-probe-profiled conformational states, as well as the rigidity and compactness of -CF3 in the target protein.

Oncogenic mutations are known to occur at non-random rates and specific hotspots, yet the structural and thermodynamic factors underlying these hotspots remain poorly understood. This study presents molecular dynamics simulations of 22 cancer hotspots across 10 oncogene families, simulated as histone-DNA complexes. Across 18 of 22 mutations, thermodynamic effects were predominantly localized to within 6 [A] of the DNA-histone interface (mean capture 104%), with van der Waals interactions driving the effect at the thermodynamic extremes, indicating that oncogenic mutations alter nucleosome binding through precise local contact changes rather than global structural rearrangements. Further, analysis revealed a gene-family correlated pattern in thermodynamic stability. RAS family mutations showed a consistent trend toward nucleosome stabilization (mean {Delta}{Delta}G = -4.82 kcal mol-1), while kinase domain mutations trended toward destabilization (mean {Delta}{Delta}G = +40.12 kcal mol-1), a difference reaching statistical significance in this exploratory analysis (Mann-Whitney U, p = 0.008). This interface-specific mechanism, combined with the gene-family-correlated thermodynamic pattern, provides a biophysical framework for understanding nucleosome-level contributions to cancer hotspot mutation biology.

Euglena longa is a heterotrophic model euglenid, of significant research interest as a unicellular eukaryote. To provide broad insight into the localisation of its proteome we used subcellular spatial proteomics (SSP), based on a label-free version of LOPIT-DC (Localisation of Organelle Proteins by Isotope Tagging after Differential Centrifugation), in which organelles are separated by differential centrifugation and each protein is assigned to a compartment from its abundance profile across the fractions. We detected 5,909 proteins, using protein abundance profiles across fractions to resolve 16 different subcellular compartments, such as the plastid, pellicle, and peroxisomes, and assigned high-confidence localisations to 2,955 of them. This study serves as a proteomic dataset to resolve protein compartmentalisation for this euglenid, enabling evolutionary, biological and functional insight and analysis. We are drafting manuscripts concerning the metabolic nature of the plastid, the evolutionary origins of the pellicle, and the functions of the peroxisome. Investigators interested in using this data prior to publication are encouraged to contact us to avoid overlap and to facilitate coordinated and collaborative use of this resource.

In our past study on human cilia from primary ciliary dyskinesia (PCD) patients using cryo-electron tomography and mass spectrometry, defects in the DNAH5 gene were proven to cause the loss of the entire outer dynein arm along with other proteins. This phenomenon was not observed in the unicellular green alga Chlamydomonas. In this study, we examined the loss of ciliary components caused by defects in various genes, based on proteomic and structural analyses of motile cilia from multiple PCD patients. The absence of single outer arm dynein genes causes a loss or significant decrease in other inner and outer arm dyneins. Interestingly, some patients with outer dynein defects show decreased intraflagellar transport (IFT) proteins, suggesting an influence of cargo components or assembly on transport. These phenomena were not observed in Chlamydomonas, suggesting a more complex mechanism of ciliogenesis in humans despite the high similarity in the final 3D architecture of the 9+2 axoneme. Defects in certain central pair proteins unexpectedly reduced components on the doublets as well. Our results demonstrate protein-protein interactions during stages of human motile ciliogenesis that are more complex than those in unicellular organisms.

Human Ca2+/calmodulin (CaM)-dependent protein kinase II alpha (CaMKII) is a highly abundant and conserved protein in the brain. It is a multifunctional serine/threonine kinase that phosphorylates cellular substrates across subcellular compartments, activating signaling pathways essential for long-term memory storage. However, a detailed understanding of protein structure and enzyme-substrate interaction mechanisms require a method to produce active, stable, and pure proteins at scale. Since phosphorylation can alter catalytic activity and substrate affinity, examining the non-phosphorylated state enables an unbiased characterization of basal enzymatic function. Therefore, producing the non-phosphorylated form of the kinase is important for biochemical and structural studies aimed at understanding kinase activation mechanisms. We describe a detailed protocol for the expression and high-yield purification of full-length, non-phosphorylated human CaMKII in bacterial cells. CaMKII was co-expressed with {lambda}-phosphatase to maintain it in a non-phosphorylated state. Protein expression was induced by autoinduction at low temperature to enhance the stability and solubility of CaMKII. Recombinant CaMKII was purified by Ni-NTA affinity chromatography, followed by gel filtration chromatography. Although CaMKII has been considered difficult to purify in its non-phosphorylated form, we obtained it in high yield, and the protein was predominantly non-phosphorylated. Our results show that combining co-expression with autoinduction improves both the solubility and yield of the homogeneous non-phosphorylated form of CaMKII, which could be suitable for subsequent biochemical and structural studies.

Extracellular vesicles (EVs) are nanosized, lipid bilayer-enclosed particles that mediate intercellular communication and hold considerable promise as biomarkers and therapeutic agents. However, their functional heterogeneity remains incompletely understood, largely due to the absence of reliable methods for isolating specific EV subpopulations. Here, we established a Nanoscale Flow sorting (NanoFACS) workflow and evaluated its potential for EV detection, sorting, and molecular characterization. We assessed critical parameters, including fluorescent staining and labeling approaches, background signal discrimination, dilution and flow-rate effects, and scatter triggering. We achieved high-purity sorting of CFSE-stained and antibody-labeled EVs, preserving vesicle integrity and sufficient RNA yields for small RNA sequencing of sorted populations. We identified EV subtype-specific small RNAs despite partial cargo overlap in single-marker sorted populations. Our results establish a reproducible antibody-based NanoFACS workflow that defines key operational parameters for reliable EV detection and sorting, with potential adaptability across conventional and spectral flow cytometry platforms. By reducing background and non-EV-associated events, the workflow enables the isolation of phenotypically defined EV subpopulations suitable for downstream molecular analysis, providing a practical approach to investigate EV heterogeneity.

Nanoplastic interactions with cell membranes can influence particle uptake, translocation, and biological effects, making their molecular description important for assessing the environmental and biological consequences of plastic pollution. However, molecular simulations generally represent nanoplastics using monodisperse polymer chains, whereas nanoplastics formed by environmental degradation are polydisperse. Here, we construct 4 nm polyethylene terephthalate (PET), polyethylene (PE), and polystyrene (PS) nanoparticles with lognormal molecular-weight distributions reported for degraded polymers, and with melt-like chain packing and entanglement. Using Martini coarse-grained molecular dynamics, we simulate their interactions for 10 microseconds with three lipid bilayers of increasing complexity. PE and PS insert rapidly and remain largely intact, whereas PET remains surface-associated and progressively releases individual chains. Fluid, polyunsaturated lipid environments undergo greater deformation and promote substantially more chain release, up to 56% of the chains for PET and 43% for PS. Release begins with the shortest chains and progressively extends to longer chains, while semicrystalline PE remains largely intact. In contrast, a size-matched monodisperse PS particle releases no chains and minimally perturbs the membrane. These results identify the molecular-weight distribution as a key determinant of nanoplastic-membrane interactions, and indicate that uniform model particles underestimate the release of polymer chains into biological membranes.

Structural model building and refinement need many decisions between different programs. These decisions often depend on an experienced researcher who knows how the programs work together and how to read experimental data. We used Claude Code with Claude Opus 5.5 to supervise this work for 24 cryo-EM and X-ray cases. The work finished in about 36 hours on one workstation. The agent prepared the inputs, ran established programs, read the validation results, compared candidate models and chose the next step. For 15 deposited cryo-EM models with weak starting geometry, the median MolProbity score improved from 2.62 to 1.76. For three crystal structures, the agent models came close to the published Rfree values. Chemical information such as ligands, ions and modifications must be provided by the researcher, in the same way as the sequence. With this information, the agent can place and check these components. These results show that agent supervision is acceptable when every decision can be inspected. The researcher makes the final scientific judgment. We suggest practical requirements for this kind of use. Future AI agents may help researchers obtain high quality and reliable structures of biological macromolecules.

Interleukin-15 (IL-15) trans-presentation (TP) by B cells is an important mechanism of T-cell activation; however, the spatial organisation of interleukin-15 receptor (IL-15R) relative to major histocompatibility complex (MHC) molecules on B-cell microvilli remains poorly understood. As microvilli protrude from the B-cell surface and may serve as sites of initial B cell-T-cell contact, the distribution of IL-15R and MHC molecules within these structures may be important during the earliest stages of T-cell recognition and activation. Here, we investigated the spatial association and molecular proximity of IL-15R with MHC class I and class II molecules on B-cell microvilli before immunological synapse formation, using confocal microscopy, stimulated emission depletion (STED) microscopy, stochastic optical reconstruction microscopy (STORM), and fluorescence lifetime imaging microscopy-based Forster resonance energy transfer (FLIM-FRET). Both MHC class I and class II molecules showed significant spatial association with IL-15R; however, the extent of colocalisation decreased as spatial resolution increased. STED microscopy revealed significant colocalisation between IL-15R and MHC class I, whereas STORM did not detect this association. In contrast, IL-15R and MHC class II remained significantly colocalised at both resolutions. FLIM-FRET further demonstrated molecular proximity between IL-15R and both MHC class I and class II molecules, with higher FRET efficiency observed for MHC class II. Collectively, these findings indicate that IL-15R is spatially organised in proximity to both MHC class I and class II molecules on B-cell microvilli before immunological synapse formation. This arrangement at potential sites of initial B-cell-T-cell contact may facilitate the coordination of IL-15 trans-presentation and antigen presentation during the earliest stages of B-cell-T-cell interactions.

Atomistic simulations of intrinsically disordered proteins (IDPs) are notoriously sensitive to force field inaccuracies, either regarding the protein or the water model, yielding inaccurate observables such as compactness, secondary structure propensities, or kinetics. The currently most widely used IDP force fields are Amber99sb-disp (A99disp) and Charmm36m (C36m). A99disp includes a new water model and thus optimized both, the protein and the water interactions. In contrast, C36m used the Tip3p water model and optimized only protein interactions. In many cases, C36m+Tip3p underestimates radii of gyration compared to FRET or SAXS experiments. Such overly compact structural ensembles are believed to arise from an imbalance between protein-protein, protein-water, and water-water interactions, which might be due to Tip3p inaccuracies. Here, we aim at re-balancing these interactions by combining C36m with the Optimal Point Charge (OPC) water model. Recently, this C36m+OPC combination showed improved accuracy for the disordered domain of the measles virus nucleoprotein. Here we present a systematic assessment, comparing C36m+OPC, C36m+Tip3p, C36m+Tip4p, C22*, A03ws, A99sb-ws, and A99disp for five IDPs, as well as a subset of those for five globular proteins, a set of disordered AGQ-repeat peptides, and the fast folding miniprotein CLN025. We compared extensive MD simulations (> 8.5 ms) with SAXS, NMR, circular dichroism, photo-induced electron transfer (PET), T-jump infrared spectroscopy, and X-ray crystallography measurements. We found that combining C36m with OPC improved accuracy over C36m+Tip3p for IDP ensembles without compromising its accuracy for globular proteins. While also the kinetics of the AGQ-peptides were more accurate for C36m+OPC, those of CLN025 folding were less accurate. Overall, C36m+OPC showed similar accuracy as A99sb-ws and A99disp, the latter is currently considered among the most accurate protein force fields.

Antifreeze proteins and peptides adhere to ice crystals and inhibit their growth, resulting in freezing point depression. AFP6 is an abundant antifreeze peptide in the blood serum of winter flounder (Pseudopleuronectes americanus). This 37-residue peptide forms an amphiphilic -helix with a hydrophobic ice-binding surface. In this study, prediction of the peptide structure suggested the formation of an antiparallel dimer burying the hydrophobic surfaces of both subunits. The oligomerisation of AFP6 was therefore examined in aqueous solutions using size-exclusion chromatography and chemical crosslinking along with ice-binding evaluation. Crosslinking of non-amidated AFP6, followed by denaturing electrophoresis, revealed bands consistent with monomers and dimers, whereas results for the amidated form and a recombinant mutant suggested additional larger forms. Bands consistent with the presence of monomers and dimers were also predominant following non-amidated AFP6 crosslinking under a variety of solution conditions, with smaller proportions of peptide assembled up to tetrameric sizes; however, the larger oligomers appeared to be absent from AFP6 treated with sodium bicarbonate. Isolated oligomers of non-amidated AFP6 produced ice crystal morphologies and thermal hysteresis comparable to those for the unfractionated peptide suggesting similar interaction with ice. Oligomerisation of AFP6 may allow greater solubility of this amphiphilic -helical peptide in aqueous solutions.

Tangential flow filtration (TFF) is an effective method for scalable EV isolation. Here, we describe and validate a reference protocol for TriPort-TFF, a laboratory-assembled TFF system featuring a three- port reservoir design enabling sample circulation, retentate return, and diafiltration. Constructed from individually sourced, commercially available components with detailed assembly and operating parameters, TriPort-TFF can be built at a lower cost than commercially integrated TFF platforms. We evaluated the system using conditioned media (CM) from human MDM cultured with EV-depleted human serum, human cerebral organoids in serum-free medium, and mouse BMDM cultured with EV-depleted foetal bovine serum (FBS) unstimulated or LPS-stimulated. The system reproducibly isolated EVs from all three sources, and biological replicates of BMDM showed distinct proteomic clusters between LPS-stimulated and unstimulated conditions. The requirement for subsequent size-exclusion chromatography (SEC) after TFF depended on the sample serum-presence. In serum-containing CM, TFF alone co-concentrated serum proteins alongside EVs, subsequent SEC improved the particle-to- protein ratio, reduced inter-replicate variability and improved proteomic depth. TFF alone in serum- free organoid CM already yielded higher particle-to-protein ratios than TFF-only preparations from serum-containing CM, with SEC further improving this ratio. TriPort-TFF offers an accessible system for EV-isolation that can be combined with SEC for diverse sample types.